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Abstract. The strawberry poison frog (Oophaga pumilio), occurring from Panama to northern Nicaragua, is characterized 
by an extraordinary variation of aposematic colouration and pattern in some Panamanian populations and has become a 
model system for studying evolutionary and ecological processes. Phylogeographic analyses of this species so far included 
populations from Panama and Costa Rica only. Based on newly collected samples from representative sites across Nicara¬ 
gua, we here extend an existing data set of mitochondrial DNA sequences of the cytochrome b gene to analyse the phylo¬ 
geographic history of O. pumilio. The 62 sequenced samples from seven localities are more closely related to the northern 
clade of O. pumilio, four haplotypes are shared amongst populations from Costa Rica and southern Nicaragua, and one 
haplotype is shared amongst Costa Rica and northern Nicaragua. An analysis of colour patterns reveals slight differences 
within and amongst Nicaraguan populations, but all of them conform to the general pattern found also in Costa Rica, with 
a reddish dorsum and blue coloured limbs. Phylogeographic analysis suggests that the Nicaraguan populations originated 
by a Late Pleistocene range expansion northwards, and reconstructs the origin of the species in southeastern Costa Rica. 

Key words. Amphibia, Anura, Dendrobatidae, Oophaga pumilio, Nicaragua, phylogeography, colouration. 


Introduction 

The Central American strawberry poison frog Oophaga 
pumilio has become an important model system for stud¬ 
ies in ecology and evolutionary biology due to its complex 
territorial and courtship behaviour (e.g., Donnelly 1989, 
Prohl & Hodl 1999, Prohl & Ostrowski 2011, Sapo- 
rito et al. 2007, Stynoski 2009), and its intriguing var¬ 
iation in aposematic colour pattern (e.g., Daly & Myers 
1967, Summers et al. 1997, Wang & Shaffer 2008, Batista 
& Kohler 2008, Maan & Cummings 2009). Populations 
of this species from the Panamanian Bocas del Toro archi¬ 
pelago are extremely polymorphic in colouration and pat¬ 
tern, amongst and in some cases also within populations 
(Daly & Myers 1967). In contrast, most mainland popula¬ 
tions are relatively uniformly coloured, with a mainly red- 
orange dorsum and at least partly blue limbs (Hagemann 
& Prohl 2007). 

The phylogeny and phylogeography of O. pumilio are 
convoluted as well. The genus Oophaga contains 9 spe¬ 
cies (O. arborea, O. granulifera, O. histrionica, O. lehmanni, 
O. occultator, O. pumilio, O. speciosa, O. sylvatica, and O. vi- 


centei (Myers et al. 1984, Grant et al. 2006)), and ranges 
from Colombia to Nicaragua. Molecular phylogenies sug¬ 
gest that O. pumilio does not form a monophyletic group 
from a mitochondrial perspective (Hagemann & Prohl 
2007, Hauswaldt et al. 2011). The mitochondrial trees 
typically separate two main lineages within O. pumilio, 
with uncorrected pairwise divergences among cytochrome 
b haplotypes of up to 7%. Furthermore, congeneric spe¬ 
cies (O. speciosa, O. arborea, O. vicentei) are phylogenet- 
ically grouped with one of these lineages (Hagemann & 
Prohl 2007, Wang & Shaffer 2008, Hauswaldt et al. 
2011). These two lineages are roughly distributed in either 
the northern or the southern part of the range, and micro¬ 
satellite data (Hauswaldt et al. 2011) also give a weak sig¬ 
nal of a northern (northeastern Costa Rica) and a south¬ 
ern (southeastern Costa Rica and Panama) cluster of pop¬ 
ulations. However, mitochondrial haplotypes belonging 
to both groups are also found co-occurring across a wide 
zone of Panama and Costa Rica. Altogether, this complex 
phylogeographic pattern suggests an intricate series of 
various past dispersal and vicariance processes, with allo- 
patric divergence and subsequent admixture, and possi- 
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ble introgression among O. pumilio and other species of 
Oophaga (Hagemann & Prohl 2007, Wang & Shaffer 
2008, Hauswaldt et al. 2011). 

All recent studies on Oophaga pumilio focused on pop¬ 
ulations from Costa Rica and Panama, omitting popula¬ 
tions from the largest portion of the distribution range 
of the species which extends across much of Nicaragua 
(Galindo Uribe et al. 2013). Only a few individuals from 
southernmost Nicaragua (Rio San Juan area) were includ¬ 
ed in general phylogenetic assessments (Grant et al. 2006, 
Santos & Canatella 2011). At present, knowledge about 
the Nicaraguan populations is scarce. As far as is known, 
they have a monomorphic reddish and blue colour pattern 
(Cope 1874, Noble 1918, Villa 1972, Savage 2002). Giv¬ 
en this lack of data, several questions about the intraspe¬ 
cific variation of O. pumilio have remained unanswered: 
(i) Do all Nicaraguan populations belong to the “northern” 
mitochondrial lineage, or are additional, strongly differ¬ 
ent haplotypes found in this part of the range? (ii) Did the 
Nicaraguan populations originate from a relatively recent 
northward range expansion? (iii) Do the Nicaraguan pop¬ 
ulations differ in colour pattern from the populations in 
northern Costa Rica? 

In this study, we aim at addressing these questions by 
studying newly collected samples of O. pumilio from seven 
localities across its range in Nicaragua. We analyse the ge¬ 
netic differentiation and phylogeography of these popula¬ 
tions based on DNA sequences of the mitochondrial cyto¬ 
chrome b gene and describe their colouration in greater 
detail. 


Material and methods 

Background: Previous knowledge on Oophaga pumilio 
from Nicaragua dates back to the early 1870s (Belt 1874, 
Cope 1874). Cope (1874) described Dendrobates ignitus 
(currently a synonym of O. pumilio) based on specimens 
collected in April 1873 at “Machuco (=Machuca) on the Rio 
San Juan, Departamento Rio San Juan, Nicaragua” (Savage 
1973). Belt (1874) recorded this species as being very abun¬ 
dant in the damp forests around Santo Domingo (Chon- 
tales Department). He experimented with feeding a young 
duck with one O. pumilio, who instead of swallowing it, 
instantly threw it up (Belt 1874). The only other contribu¬ 
tion on toxins of O. pumilio from Nicaraguan populations 
is that of Mebs et al. (2008), who studied alkaloids in pop¬ 
ulations from Cerro Musun in the central part, and along 
the Rio San Juan and Rio Indio in the southeastern portion 
of the country. The Musun population exhibited a largely 
different alkaloid profile compared with that of Panama¬ 
nian specimens, the latter sharing the common alkaloids 
with those from Rio San Juan/Rio Indio (Mebs et al. 2008). 
However, these profiles are of no direct phylogenetic rel¬ 
evance since alkaloids in dendrobatids are of dietary origin 
(Mebs et al. 2008). Only a few other data are available on 
the ecology and biology of the species in Nicaragua: Noble 
(1918) studied the stomach contents of several O. pumilio 
collected around Cuckra (today in the Atlantico Sur De¬ 


partment), finding mostly small red ants, a single pill bug, 
and a spider, and Wong et al. (2009) studied heterospecific 
acoustic interference effects on the calling of adult males 
O. pumilio at Bartola, Rio San Juan Department. 

Most of the historical distribution records of this species 
from Nicaragua still are represented by extant populations, 
although these are often reduced in size. For example, the 
populations mentioned by Belt (1874) today persist in a 
small forest patch close to the Reserva Privada Las Brumas, 
but are threatened by deforestation. Also the populations 
referred to by Cope (1874) are still found throughout the 
southeastern portion of the country. Noble (1918) men¬ 
tioned populations that are now located in a small forest 
patch in the Reserva Natural Kahka Creek that was spared 
the devastating effects of hurricane Felix in 2007. On the 
other hand, the historical populations around the city of 
Matagalpa (Matagalpa Department) have not been con¬ 
firmed since the mid 1960s, and there are other uncertain¬ 
ties regarding the historical distribution of O. pumilio in 
Nicaragua. As an example, Gaige et al. (1937) did not re¬ 
port O. pumilio from an extensive three-month herpeto- 
logical collection carried out in 1935 in the area around the 
Rio Escondido and its tributaries in central-eastern Nica¬ 
ragua (today in the Atlantico Sur Department), an appar¬ 
ently suitable area, while the species is known to occur to 
the north, south, east and west of this region. Additionally, 
there have been several questionable records of O. pumilio 
in Nicaragua from the surroundings of Siuna in the Atlan¬ 
tico Norte Department (Villa 1972) as well as from Mana¬ 
gua and Chinandega Departments. It is possible that the 
voucher specimens for the Atlantico Norte records, as well 
as several other sites shown in the distribution map of Vil¬ 
la (1972), were lost during the catastrophic earthquake in 
Managua in 1972, which destroyed all of the alcoholic her- 
petological specimens housed in the Museo Nacional de 
Nicaragua (Villa 1981). 

Currently, the distribution of O. pumilio in Nicaragua 
is considerable compared to its entire range. In Nicaragua, 
O. pumilio occurs across a continuous area in the south¬ 
eastern corner of the country (Departamento Rio San 
Juan; Kohler 2001, Sunyer et al. 2009) and other remain¬ 
ing lowland forest patches of the Departamento Atlantico 
Sur. Furthermore, the species is found in a few isolated and 
threatened forest localities at mid-altitudes in the Depar- 
tamentos Chontales, Boaco, and Jinotega (Kohler 2001). 
In Nicaragua, this species is found in association with run¬ 
ning waters and forests with different degrees of anthro¬ 
pogenic alteration. Individuals sometimes can be found 
in high densities, even in grasslands where no substrates 
for refuge and oviposition are obvious. However, it is un¬ 
known how these populations respond to increased habitat 
degradation and transformation. 

Within Nicaragua, O. pumilio is considered a species of 
medium vulnerability and known to occur in four forest for¬ 
mations and several bioclimatic regions, with a widespread 
distribution in Lowland Wet and a peripheral distribu¬ 
tion in Lowland Moist, Premontane Wet, and Premontane 
Moist formations, and can be found from sea level to 960 m 
a.s.l. (Sunyer & Kohler 2010). According to the IUCN cri- 
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teria, this species is regarded as a species of Least Concern, 
although populations are decreasing (Soli's et al. 2013), and 
according to the Convention on the International Trade of 
Endangered Species (CITES), it has been considered an Ap¬ 
pendix II species since 1987 (CITES 2013). Oophaga pumilio 
is one of the most profitable amphibian species for the pet 
trade in Nicaragua, and between 1991 and 1996, 95% of the 
exports of this species were shipped to the USA (Soli's et 
al. 2013). Even though Cerro Kilambe is not included in the 
IUCN distribution map of the species, this locality was well 
known to the pet traders in the early 1990s who offered lo¬ 
cals one US dollar per live specimen (pers. comm, by local 
people to J. Sunyer): a large amount in the aftermath of an 
impoverishing civil war. Nicaragua then fixed export quo¬ 
tas to 1,100 specimens taken from the wild in 1997, 1,200 


ranched specimens in 1998, 10,000 ranched specimens in 
1999, and 3,450 ranched and captive-bred specimens in 
2000 and 2001, respectively (CITES 2013). 

Sampling localities and procedures 

Specimens of Oophaga pumilio were collected from seven 
localities in Nicaragua (Fig. 1. Tab. 1) between May and Au¬ 
gust 2012. In the northernmost distribution area, we sam¬ 
pled two highland populations from the Reservas Naturales 
Cerro Kilambe (Jinotega region) and Cerro Musun (Mata- 
galpa region), located at higher altitudes in Nicaragua (770 
and 572 m a.s.l., respectively). The Kilambe population also 
represents an extension of the known range and currently 


85°0'0"W 84°0'0"W 83°0'0"W 82°0'0"W 



Figure 1. Location of the Nicaraguan populations of O. pumilio included in this study, and of the Costa Rican and Panamanian popula¬ 
tions previously studied by Hauswaldt et al. (2011). Numbers correspond to localities listed in Table 1. 
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Table 1. Locality numbers and names, coordinates, altitudes, and numbers of individuals of 0. pumilio sequenced for the cytochrome b 
fragment, from Nicaraguan (NIC) populations (this study), and from Costa Rican (CR) and Panamanian (P) populations previously 
studied by Hauswaldt et al. (2011). Populations and numbers of individuals (in parenthesis) assigned to the 93 haplotypes (H) found 
in 321 individuals of O. pumilio, two of Oophaga vicentei (O.v.), and one of Oophaga granulifera (O.g.). 


Locality Locality Province Coordinates [lat. Cyt b Haplotypes 

number (°N), long. (°W)| (N) 


1 

NIC: Kilambe 

Jinotega 

13.63 

-85.72 

8 

H5 (7), H6 (1) 

2 

NIC: Musun 

Matagalpa 

12.96 

-85.23 

8 

HI (7), H5 (1) 

3 

NIC: Kahka Creek 

Atlantico Sur 

12.67 

-83.72 

7 

H5 (7) 

4 

NIC: Las Brumas 

Chontales 

12.28 

-85.08 

9 

H6 (9) 

5 

NIC: Guatuzos 

Rio San Juan 

11.01 

-85.06 

10 

H9 (10) 

6 

NIC: Bartola 

Rio San Juan 

10.97 

-84.34 

10 

HI (1), H2 (2), H7 (2), H8 (5) 

7 

NIC: Rio Indio Lodge 

Rio San Juan 

10.93 

-83.73 

10 

HI (1), H2 (3), H3 (1), H4 (5) 

8 

CR: Upala 

Alajuela 

10.91 

-85.05 

7 

H9 (4), H19 (3) 

9 

CR: Cano Negro 

Alajuela 

10.87 

-84.78 

11 

H20 (5), H21 (1), H22 (2), H23 (3) 

10 

CR: Tortuguero 

Limon 

10.61 

-83.53 

9 

HI (4), H8 (2), H10 (3) 

11 

CR: La Selva 

Heredia 

10.43 

-84.00 

11 

H2 (2), H8 (1), H33 (1), H34 (2), H35 (3), H36 (1), 
H37 (1) 

12 

CR: Pueblo Nuevo 


10.32 

-83.59 

11 

H8 (2), H24 (7), H25 (1), H26 (1) 

13 

CR: Guapiles 

Limon 

10.19 

-83.82 

11 

H15 (1), H27 (2), H28 (8) 

14 

CR: Siquirres 

Limon 

10.10 

-83.52 

12 

H14 (1), H15 (1), H16 (6), H17 (2), H18 (2) 

15 

CR: Rio Reventazon 

Cartago 

10.09 

-83.56 

9 

H16 (4), H27 (1), H38 (2), H39 (1), H40 (1) 

16 

CR: Hitoy Cerere 

Limon 

9.67 

-83.09 

14 

H8 (1), H27 (1), H58 (6), H89 (4), H90 (2) 

17 

CR: Puerto Viejo 
de Talamanca 

Limon 

9.65 

-82.76 

12 

Hll (5), H12 (5), H13 (2) 

18 

CR: Bribri 

Limon 

9.65 

-82.88 

9 

H29 (3), H30 (3), H31 (1), H32 (2) 

19 

P: Colon 

Bocas del Toro 

9.39 

-82.24 

12 

H41 (8), H42 (1), H43 (1), H44 (1), H45 (1) 

20 

P: Solarte 

Bocas del Toro 

9.33 

-82.22 

12 

H42 (1), H50 (11) 

21 

P: Bastimentos 

Bocas del Toro 

9.30 

-82.14 

10 

H41 (3), H42 (1), H48 (5), H49 (1) 

22 

P: Almirante 

Bocas del Toro 

9.29 

-82.39 

9 

H50 (3), H51 (2), H52 (3), H53 (1) 

23 

P: San Cristobal 

Bocas del Toro 

9.27 

-82.29 

10 

H8 (1), H42 (2), H83 (4), H84 (3) 

24 

P: Pastores 

Bocas del Toro 

9.24 

-82.35 

10 

H78 (2), H79 (4), H80 (2), H81 (1), H82 (1) 

25 

P: Pop a 

Bocas del Toro 

9.22 

-82.13 

19 

H46 (1), H50 (1), H69 (7), H70 (1), H71 (4), H72 (1), 
H73 (1), H74 (1), H75 (1), H76 (1) 

26 

P: Tierra Oscura 

Bocas del Toro 

9.18 

-82.26 

8 

H54 (5), H56 (2), H77 (1) 

27 

P: Punta Alegre 

Bocas del Toro 

9.16 

-81.91 

8 

H85 (3), H86 (3), H87 (1), H88 (1) 

28 

P: Cayo de Agua 

Bocas del Toro 

9.16 

-82.05 

12 

H50 (7), H59 (1), H60 (1), H61 (1), H62 (1), H63 (1) 

29 

P: Cauchero 

Bocas del Toro 

9.16 

-82.25 

12 

H53 (3), H54 (2), H55 (1), H56 (4), H57 (1), H58 (1) 

30 

P: Loma Partida 

Bocas del Toro 

9.14 

-82.17 

10 

H46 (9), H47 (1) 

31 

P: Escudo de Veraguas Escudo de Veraguas 

9.10 

-81.55 

11 

H64 (7), H65 (1), H66 (1), H67 (1), H68 (1) 






321 



O.v.l 

Veraguas 

8.51 

-81.08 

1 

H91 (1) 


O.v.2 

Code 

8.63 

-80.58 

1 

H92 (1) 


o.g. 




1 

H93 (1) 


is the northernmost known locality of the species. In cen¬ 
tral Nicaragua, we sampled frogs from the private Reserva 
Natural Las Brumas at 562 m a.s.l., and from Kahka Creek, 
a lowland (24 m a.s.l.) population near Laguna de Perlas on 
the Atlantic coast. In the southern part of Nicaragua, we 
sampled three populations close to the border with Costa 
Rica: the Refugio de Vida Silvestre Los Guatuzos close to 
the Papaturro River, the Reserva Silvestre Privada Refugio 


Bartola, and the Rio Indio Lodge near the Atlantic coast, at 
54, 74 and 14 m a.s.l., respectively. Geographic coordinates 
were obtained with a Garmin Oregon 550 GPS. 

Ventral and dorsal photographs were taken with a dig¬ 
ital camera of ten specimens per site. To facilitate body 
measurements, every frog was placed on a millimetre grid. 
Colour descriptions were made considering mainly the ter¬ 
minology used by Savage (1968, 2002). Snout-vent length 
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(SVL) measurements were taken based on the photographs 
using ImageJ r.46r (Abramoff et al. 2004). Tissue samples 
were taken by toe-clipping from up to ten specimens per 
locality and preserved in 96% ethanol. 

DNA sequencing 

Genomic DNA was extracted from toe clippings using the 
QIAGEN DNeasy Blood and Tissue kit. Samples and ex¬ 
tractions were stored at -20°C. We amplified a fragment of 


the mitochondrial cytochrome b gene using the primers 
MVZ r5 (Moritz et al. 1992), and CytbDenr-H (Santos 
& Cannatella 2011). PCR conditions consisted of an ini¬ 
tial 2 min. at 95°C, 35 cycles of 30 sec. each at 95°C, 40 sec. 
at 42°C, and 60 sec. at 72°C, and a final extension step of 
6 min. at 72°C. Amplified fragments were sequenced with 
the forward primer MVZ15, using dye-terminator chem¬ 
istry. DNA sequences (784 bp) were subsequently aligned 
with those from Hauswaldt et al. (2011). All newly de¬ 
termined sequences were deposited in GenBank (accession 
numbers KF645289-KF645350). 



Figure 2. Oophaga pumilio colour patterns observed at Nicaraguan sites. 1 - Kilambe; 2 - Musun; 3 - Kahka Creek (Photo taken from 
www.monarchzman.deviantart.com); 4 - Las Brumas; 5 - Guatuzos; 6 - Bartola; 7 - Rio Indio Lodge. 
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Table 2. Dorsal and ventral colours and patterns described for Oophaga pumilio populations in Nicaragua. Colour pattern descriptions 
and snout-vent length (SVL) measurements (± standard deviation) were taken for 10 individuals. 


Population 

Dorsal body Dorsal limb colour Dorsal body 
colour pattern 

Dorsal leg 
pattern 

Ventral body 
colour 

Ventral 

limb 

colour 

Ventral body 
and limb 
pattern 

SVL 

Kilambe 

Red 

Blue with or with¬ 
out red markings 
on one thigh 

None or very 
few scattered 
small spots or 
wavy flecks 

Reticulated 

markings 

Mainly red 
with light blue 

Blue 

Diffused flecks 

22.4±0.6 

Musun 

Red-orange 

Blue, with or with¬ 
out red markings 
on one thigh 

Small spots or 
wavy flecks 

Reticulated 
markings or 
diffused spots 

Mainly orange- 
red with light 
blue 

Blue 

Diffused small 
spots or flecks 

22.7±1.1 

Kahka 

Creek 

Red and 
blue near 
cloaca 

Blue 

Small spots or 
wavy flecks 

Diffused reticu¬ 
lated markings or 
small spots 

Mainly light 
blue with red 
in the flanks 

Blue 

Diffused flecks 

22.2±0.8 

Las Brumas 

Red or 
orange red 

Blue, with or with¬ 
out red markings 
on thighs 

Small spots or 
wavy flecks 

Mainly well-de¬ 
fined reticulated 
markings 

Mainly orange- 
red with light 
blue 

Blue 

Diffused flecks 

20.9±0.8 

Guatuzos 

Red 

Blue 

None or very 
few scattered 
small spots or 
wavy flecks 

Small spots 

Mainly orange- 
red with light 
blue near cloaca 
and chin 

Blue 

None or a few 
diffused flecks 

23.4±0.7 

Bartola 

Dull red 

Grey-blue with or 
without red mark¬ 
ings on thighs 

Very few 
scattered 
small spots 

Diffused spots 

Orange-red 
with light blue 
or red with light 
blue near cloaca 

Light 

blue 

Diffused flecks 

19.7±0.8 

Rio Indio 
Lodge 

Red-orange 
or orange- 
red 

Red-orange, with 
or without black 
or blue at tips 

None or small 
spots or wavy 
flecks 

Small spots or 
stripes 

Orange red 

Red, 
black, 
and blue 

Diffused flecks 

21.8±1.1 


Phylogeography 

The Bayesian Information Criterion (BIC) incorporat¬ 
ed in jModeltest vs. 2.1.1 (Darriba et al. 2012) was used 
to determine the best fitting model of nucleotide substitu¬ 
tion for two datasets of the cytochrome b sequences. The 
HKY + T model was selected for dataset 1 composed of all 
haplotypes of O. pumilio as well as of the outgroup species 
(O. granulifera and O. vicentei ), and the TrN + T model was 
selected for the dataset 2 containing sequences of all indi¬ 
viduals of O. pumilio. 

Dataset 1 was used to reconstruct a phylogenetic tree 
of haplotypes using Bayesian inference as implemented in 
Mr. Bayes vs. 3.1.2. (Ronquist & Huelsenbeck 2003) us¬ 
ing an unpartitioned approach. Using the default settings 
of the software, we ran two chains for 10 million genera¬ 
tions each, sampling every 100 th tree, and checked for con¬ 
vergence. The first 25% of the sampled trees were discarded 
as burn-in. Bayesian posterior probabilities are based on a 
majority-rule consensus tree. 

Dataset 2 was used with the software Beast 1.7. (Drum¬ 
mond et al. 2012) to simultaneously estimate tree topo¬ 
logy and geographic locality of nodes, using the geograph¬ 
ic coordinates of the samples as traits and combined this 
with the phylogeographic method by Lemey et al. (2010). 
Specifically, we used the Relaxed Random Walk model 
to estimate the geographic origin of the O. pumilio clade 


and possible dispersal routes throughout Panama, Costa 
Rica and Nicaragua. We set the coalescent prior to con¬ 
stant size and used an uncorrelated relaxed clock model 
(Drummond et al. 2006), the TrN + T model of substi¬ 
tution, codon partitioning (1+2, 3), and set the mutation 
rate to 0.01. We ran a MCMC chain 3 x 10 8 long, sampling 
every 30.000 th . We ensured that effective sample sizes 
(ESS) were higher than 200 for all parameters with Tracer 
1.5. The first 20% of the trees were discarded as burn-ins. 
The resulting tree was summarized with Tree Annotator 
1.7.x. To generate the KML file that could be plotted with 
Google Earth (www.googleearth.com), we used the soft¬ 
ware SPREAD (Bielejec et al. 2011) and a posterior prob¬ 
ability setting of 0.05. 

Results 

Morphology and colour pattern 

No substantial variation in colour pattern and body size 
was detected among Nicaraguan O. pumilio populations. 
All share the “blue jeans” dorsal colour pattern seen in Cos¬ 
ta Rican frogs, with blue or black limbs (Figs 2-3, Tab. 2). 
The colour of the body differs slightly amongst populations 
and ranges from dull red to red-orange. The colour of the 
limbs varies from blue to grey-blue or black. The presence 
and amount of black on the dorsum can differ as well and 
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take on the shape of small spots, wavy lines, or reticulated 
markings. The venter is orange or red, with or without light 
blue (Fig. 3). Specimens from Kilambe and Guatuzos have 
a particularly bright red body colour and blue limbs with 
reduced black pattern on the dorsum, and those from Gua¬ 
tuzos have very small light blue markings on the venter. 
Individuals from Musun have a red-orange dorsum with 
a black pattern. Specimens from Kahka Creek have a red 
dorsum with some blue on the cloaca and limbs, and a re¬ 
duced amount of red on the venter and flanks. Specimens 
from Bartola have a dull red-coloured dorsum with light 
blue limbs, and are mainly red on the ventral side. At Las 
Brumas, the frogs are either red or red-orange with blue 
legs, and always exhibit a black dorsal pattern. At Rio Indio 
Lodge, two main morphs can be distinguished: red-orange 
and blue-legged frogs with a black dorsal pattern, and or¬ 
ange-red frogs with or without black legs and a very re¬ 


duced or even absent dorsal pattern (Fig. 2, Tab. 1). Popula¬ 
tions varied in mean body size between 19.7 and 22.7 mm, 
with the largest frogs occurring in the Rio San Juan popu¬ 
lations from Guatuzos, and the smallest frogs being found 
in Bartola (Tab. 2). 


Phylogeography 

A total of 62 cytochrome b sequences were obtained 
from Nicaraguan O. pumilio and compared with previ¬ 
ously determined sequences from Costa Rica and Pana¬ 
ma (Hauswaldt et al. 2011) to yield a total alignment of 
sequences from 321 individuals. Because the Costa Rica 
and Panama sequences were shorter than the newly deter¬ 
mined ones from Nicaragua, the alignment used for analy¬ 
sis was only 558 bp in length. 88 haplotypes could be distin- 
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Figure 3. Representation of the terminology used for describing the dorsal and ventral colours and patterns of Oophaga pumilio 
populations from Nicaragua. 
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Figure 4. MrBayes tree 
based on an unparti¬ 
tioned dataset of 529 bp 
cyt-b sequences (all 
haplotypes are of equal 
length), using the 
HKY+G model and 
O. granulifera as an out¬ 
group. Red branches and 
numbers indicate haplo¬ 
types found in Nicara¬ 
gua. Grey scales in boxes 
represent the relative fre¬ 
quency of individuals as¬ 
signed to each haplotype, 
from the region of North 
Nicaragua (NN), Central 
Nicaragua (CN), South 
Nicaragua (SN), Costa 
Rica (CR), and Panama 
(PA). 
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guished (Fig. 4). The phylogenetic analysis (Fig. 4) agrees 
with previous findings in separating the haplotypes into a 
northern and a southern clade, and in placing sequences of 
O. vicentei with the northern clade of O. pumilio. A total of 
nine haplotypes (numbers as in Fig. 4), all grouped in the 
northern clade, were found among the Nicaraguan popula¬ 
tions. Five of these (H3—H7) had not been observed before 
and were exclusive to Nicaraguan populations. The other 
four haplotypes (Hi, H2, H8, H9) were shared among Nic¬ 
araguan and Costa Rican populations. Haplotype sharing 
with Costa Rican populations was not limited to the popu¬ 
lations from Rio San Juan (geographically closest to Costa 
Rica), but was observed in northern Nicaraguan popula¬ 
tions as well: e.g., haplotype Hi, was found at Musun, Bar- 
tola, Rio Indio Lodge, and in the Costa Rican population 
ofTortuguero (Tab. 1). 

The phylogeographic analyses strongly support a cen¬ 
tre of origin of O. pumilio in southeastern Costa Rica from 
where the species spread north- and southwards. Disper¬ 
sal to the island of Escudo de Veraguas likely co-occurred 
with the spread into Nicaragua (Fig. 5). The origin of the 
populations in Nicaragua is reconstructed as a recent Pleis¬ 
tocene dispersal at approximately 200,000 years before 
present (Fig. 5). 


Discussion 

Understanding the phylogeographic history of O. pumilio 
is crucial for making inferences on the evolution of the ex¬ 
traordinary polychromatism of the species. However, giv¬ 
en the lack of molecular data from Nicaragua, i.e., from 
more than half of the species’ range, it used to be impos¬ 
sible to reliably test hypotheses on the geographic origin 
and directions of range expansion of this species. The data 
presented herein partly fill this geographical sampling 
gap, and confirm that the Nicaraguan populations are not 
strongly differentiated, but rather the result of a recent 
range expansion from Costa Rica. Our explicit phylogeo¬ 
graphic analyses place the origin of the species in southern 
Costa Rica, and support a radial expansion of this species 
northwards into Nicaragua and southwards towards the 
Bocas del Toro archipelago during the Late Pleistocene. 
However, an important limitation of this analysis is its sole 
reliance on mitochondrial DNA. If either the northern or 
southern lineage of O. pumilio were the result of an ancient 
hybridisation and introgression process with another spe¬ 
cies of Oophaga, then some of the detailed reconstructions 
herein might not fully reflect the true evolutionary history 
of the species. Despite this limitation, given that the north- 



Figure 5. Phylogeographic reconstruction for Central American Oophaga pumilio populations using the Relaxed Random Walk Model. 
Green polygons and red lines indicate early events, whereas black lines and polygons more recent events. The red marking points out 
the Isla Bastimentos in Panama. MA - million years ago. 
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ern Nicaraguan populations of O. pumilio share haplotypes 
with Costa Rican populations, and that Nicaraguan popu¬ 
lations overall do not differ strongly in colouration from 
the ones in Costa Rica, we find it very plausible that Nica¬ 
raguan populations indeed originate from a recent range 
expansion northwards. 

The relationships among mitochondrial lineages and 
species of Oophaga also reveal the need for taxonomic revi¬ 
sion. Daly & Myers (1967) and Myers & Daly (1976) sug¬ 
gested recognizing O. pumilio as a single polymorphic spe¬ 
cies, and this view has either been accepted or controver¬ 
sially discussed (Summers et al. 1997, Hagemann & Prohl 
2007, Batista & Kohler 2008, Hauswaldt et al. 2011). 
The taxonomic integrity of O. pumilio as a single entity re¬ 
mains uncertain, especially relative to several other species 
in the genus distributed in Panama (O. arborea, O. speciosa, 
and O. vicentei ). Alternatively, some (or all) of these oth¬ 
er taxa could in fact represent morphs of O. pumilio and 
thus junior synonyms, or some populations considered as 
O. pumilio could in fact rather belong to one of those close¬ 
ly related species. The strawberry poison frog was original¬ 
ly described as Dendrobates pumilio Schmidt, 1857 from 
the Bocas del Toro province in Panama. Daly & Myers 
(1967) applied this name to the small-sized poison frogs 
occurring from Nicaragua to northwestern Panama, and 
Grant et al. (2006) suggested classifying this and related 
species, separate from Dendrobates, in the genus Oophaga. 
Three other names are currently considered to be junior 
synonyms of O. pumilio: Dendrobates galindoi Trapido, 
1953, described likewise from Bocas del Toro, Dendrobates 
typographus Keferstein, 1867, and Dendrobates ignitus 
Cope, 1874, with their type localities in Costa Rica and Nic¬ 
aragua, respectively (see Savage 1968, Batista & Kohler 
2008). However, given the wide co-occurrence of the two 
haplotype lineages in O. pumilio (Hauswaldt et al. 2011), 
we consider it to be rather unlikely that they represent dis¬ 
tinct and independent evolutionary lineages, and that any 
of these junior synonyms would require elevation to spe¬ 
cies level. Inversely, as mentioned above, the identities, sta¬ 
tus and distribution ranges, especially of O. arborea and 
O. speciosa, but also of O. vicentei, require detailed revision. 
A first step that has to be taken towards such a revision is 
a comprehensive mitochondrial phylogeny of all Oophaga 
species, using longer DNA sequences, in order to obtain a 
strongly supported phylogenetic tree depicting the mito¬ 
chondrial relationships among these frogs. Using phyloge- 
nomic methods such as RAD sequencing (e.g., Rubin et al. 
2012) would then allow to reconstruct the true evolution¬ 
ary relationships among species and populations, and to 
understand if, when, and where phenomena of mitochon¬ 
drial introgression might have taken place in this genus. 
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